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Recent experimental breakthrough in realizing spin-orbit (SO) coupling for cold atoms has spurred 
considerable interest in the physics of 2D SO coupled Fermi superfluids, especially topological Ma- 
jorana fermions (MFs) which were predicted to exist at zero temperature. However, it is well known 
that long-range superfluid order is destroyed in 2D by the phase fluctuation at finite temperature 
and the relevant physics is the Berezinskii-Kosterlitz-Thouless (BKT) transition. In this Letter, we 
examine finite temperature effects on SO coupled Fermi gases and show that finite temperature is 
indeed necessary for the observation of MFs. MFs are topologically protected by a quasiparticle en- 
ergy gap which is found to be much larger than the temperature. The restrictions to the parameter 
region for the observation of MFs have been obtained. 
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A new research direction in low dimensional condensed 
matter physics that attracts much recent attention is 
the study of two-dimensional (2D) topological quantum 
states of matter (e.g., fractional quantum Hall effects, 
chiral p-wave superfluids/super conductors, etc.) that 
support exotic quasiparticle excitations (named anyons) 
with Abelian or non-Abelian exchange statistics [lMs|- 
For instance, it has been shown recently that a topolog- 
ical cold atom superfluid may emerge from an ordinary 
2D 5-wave Fermi superfluid in the presence of two addi- 
tional ingredients: spin-orbit (SO) coupling and Zeeman 
field [9]. Such topological superfluids can host Majo- 
rana fermions (MFs), non-Abelian anyons which are their 
own antiparticles, and may have potential applications in 
fault-tolerant topological quantum computation 

On the experimental side, 2D degenerate s- wave 
Fermi gases have been realized using highly anisotropic 
pancake-shaped trapping potential jlCH13|. Further- 
more, SO coupling and Zeeman field for cold atoms have 
been generated in recent pioneer experiments through the 
coupling between cold atoms and lasers [I3-[l9| . It seems 
therefore that MFs are tantalizingly close to experimen- 
tal reach in degenerate Fermi gases. Inspired by the ex- 
perimental achievements, there have been extensive the- 
oretical efforts for understanding the physics of SO cou- 
pled Fermi gases, particularly the mean- field crossover 
from the Bardeen-Cooper-Schrieffer (BCS) superfluids to 
Bose-Einstein condensation (BEC) of molecules |2CH33j. 
While the mean-field theory works qualitatively well in 
3D, it may not yield relevant physics in 2D at finite tem- 
perature. For instance, differing from the mean-field pre- 
diction, there is no long-range superfluid order in 2D at 
finite temperature due to phase fluctuations At fi- 
nite temperature, the relevant physics is the BKT transi- 
tion [35|,l36[ with a characteristic temperature Tbkt, be- 
low which free vortex-antivortex (V-AV) pairs are formed 



spontaneously. When the temperature is further lowered 
below another critical temperature Tvortex, the system 
forms a square V-AV lattice [37|, [38| . Because MFs only 
live in the cores of spatially well separated vortices, it is 
crucial to study the dependence of Tbkt and Tv or tex on 
the SO coupling strength in 2D Fermi superfluids. 

In this Letter, by taking account of phase fluctuations 
in the finite temperature quantum field theory, we inves- 
tigate the finite temperature properties of SO coupled 
degenerate Fermi gases en route to clarifying some cru- 
cial issues for the experimental observation of MFs in 
such systems. Our main results are the following: 

(I) We show that, quite unexpectedly, the SO coupling 
reduces both Tv or tex and Tbkt, despite it enhances the 
superfluid order parameter A. 

(II) In the pseudogap phase (T > Tbkt), the Fermi 
gas lacks phase coherence. While in the vortex lattice 
phase (T < Tv or tex), the short distance between neigh- 
boring vortices may induce a large tunneling between vor- 
tices that destroys the Majorana zero energy states [3^ . 
Therefore MFs can be observable only in the cores of nat- 
urally present V-AV pairs [40] in the temperature region 
Tvortex < T < Tbkt, instead of intuitively expected zero 
temperature. 

(III) MF in a vortex core is protected by a quasipar- 
ticle energy gap > A 2 /2Ef (Ep is the Fermi energy), 
which is found to be much larger than the experimental 
temperature. Such a large energy gap greatly reduces the 
occupation probability of non-topological excited states 
in the vortex core due to finite temperature, and ensures 
the topological protection of MFs. 

(IV) The restrictions to the parameter region for the 
observation of MFs have been obtained. We also propose 
a simple strategy for finding experimental parameters for 
the observation of MFs. 

We consider a 2D degenerate Fermi gas in the presence 
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of a Rashba type of SO coupling and a perpendicular Zee- 
man field. In experiments, 2D degenerate Fermi gases 
can be realized using a ID deep optical lattice with a po- 
tential Vq sm 2 (27r z/a w ) along the third dimension, where 
the tunneling between different layers is suppressed com- 
pletely [10Ml3|. The Rashba SO coupling and Zeeman 
field can be realized using adiabatic motion of atoms in 
laser fields [H], [Hj • The Hamiltonian for this system can 
be written as (H = Kb = 1) 



H = H* 



(1) 



where the single atom Hamiltonian Hp = a= ^ ^(ek — 
/ i c r )^kcr^ kcr ' ^kcr i s ^ ne crea tion operator for a fermion 



atom with momentum k and spin <j, ek = /c 2 /2m, 
m is the atom mass, ^ = \i + /i, = fi — h, fi 
is the chemical potential, and h is the Zeeman field. 
The Hamiltonian for the Rashba type of SO coupling 
is H soc = aEklft ~ ik x )Cl t C^ + (k y + ik^C^C^}. 
The interaction between atoms is described by Hi = 
~9 Xlk ^-kt^k|^ k i^- k t> where the effective regularized 
interaction parameter 1/q = ^ k l/(2ek + E^) for a 2D 



Fermi gas [28|, [3l|, [32|, |4l|. In experiments, the bind- 



ing energy E}> can be controlled by tuning the s-wave 
scattering length or the barrier height Vo along the z di- 
rection. Small and large E^ correspond to the BCS and 
BEC limit, respectively [Hj. 

The finite temperature properties of the 2D Fermi 
gas are obtained using finite temperature quantum field 
theory, where the action for the Hamiltonian (pQ) is 
Sv = JodT[E Ka Ci a d T C^ + H] with = 1/T. In- 
troducing the standard Hubbard- Stratonovich transfor- 
mation with the mean field superfluid order parameter 
(j) = g J2k (Ck^C-kt) an d integrating out the fermion de- 
grees of freedom, we have the partition function Z = 
J D(\)D(\f exp(— S^ff) with the effective action S e s = 
Jo drig- 1 |0| 2 + Ck) - f TVpnG" 1 ]. Here 
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is the inverse Nambu matrix under the Nambu basis 
¥(k) = (C kT ,C H ,Cl H ,Cl kT ) T , the symbol Tr de- 
notes the trace over momentum, imaginary time, and 
the Nambu indices, k± = (k y ± ik x ) and C±M ~ ( 6k — 
/i) zb h = ("k ± h. The superfluid order parameter A is ob- 
tained from the saddle point of the effective action (i.e., 
=A = 0), which yields the gap equation 
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Here the quasiparticle energy spectrum £?±,k = 
y/E£ + a 2 k 2 + h 2 ± 2 A, A = ^/a 2 ( 2 k 2 + h 2 E 2 , A± = 
\{l ± h 2 /A), E 2 =Cl + |A| 2 , and fc 2 = k 2 x + k 2 y . 



In 2D Fermi gases, it is well known that long-range 
superfluid order can be destroyed by phase fluctuations 
of the order parameter at any finite temperature [34J . To 
study the phase fluctuations in SO coupled Fermi gases, 
we set (j) = Ae z9 following the standard procedure, where 
6 is the superfluid phase around the saddle point (deter- 
mined by Eq. (|2j)) that varies slowly in position and time 
spaces. The superfluid phase can be decoupled from the 
original Green's function through a unitary transforma- 



tion UG- 1 (9)U ] = Gq where U 
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= exp(zM<9/2), and 
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M = diag(l, 1,-1,-1). ^- ,, v 2 , gm 
^m V ) ^( T ^ CF xdyO — Icr y d x 9) is the corresponding self- 
energy, <Ji and Ti are Pauli matrices in the Nambu space. 
Gq 1 is the Green's function at <p = A. The effective ac- 
tion can be written as S^ff = -So (A) + 5fl U c(V#, d T 0) j43[, 
where S fluc (V#, d T 6) = trE n >i H G o^) n - Expanding £ 
up to the leading order (n = 2), we have 
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Xlk nk - ^ e nave checked 
that the expressions for J and Q reduce to previous 
results in different limits j4ll. E3 - [46| . 

The phase 9 can be decomposed into a static vortex 
part # v ( r ) and a time-dependent spin- wave part sw (r, r). 
As a consequence, the action of the phase fluctuation 
becomes 5fl uc = S v + 5 SW with S v = \ J d 2 r J[v6>(r)] 2 
and ^ sw = i/d 2 r{J[V# sw (r,r)] 2 + P[S r ^ sw (r, r)] 2 - 
Q[i9 r 6> sw (r, r)]} = X^k 111 ! 1 exp(-^cj k )], where cj k = 
c |k|, c = \fj~JP is the speed of the spin wave |44|. Note 
that Eq. (j3j) is exactly the same as the effective action 
for the 2D Heisenberg XY model 35, [36| with different 
J, P and Q, therefore the BKT transition temperature 
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Across the BKT temperature, there is a transition from 
the pseudogap phase (with finite pairing A but without 
phase coherence or superfludility) to phase coherent V- 
AV pairs (with both pairing and superfluidity). When 
the temperature is further lowered, free V-AV pairs form 
a tightly bounded V-AV lattice below another critical 



temperature 



37, 
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The atom density equation can be obtained from the 
total thermodynamic potential Q = TS G s, yielding [i?}, 



n = —dQ/dfi = \^ rik — f3 1 dS sw /dfi, 



(6) 
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FIG. 1: (Color online) The dependence of the BKT tempera- 
ture on SO coupling (a), Zeeman field (b), and binding energy 
(c,d). ft = in (a) and (c). aK F = 1.0E F in (b) and (d). 



where the atom density n = mE^/n, and Ep = 
h 2 Kp/2m is the Fermi energy without SO coupling and 
Zeeman field. The length unit is chosen as the inverse of 
the Fermi vector Kp 1 . 

We numerically solve Eqs. (j2j), (j4]) and (j6j) self- 
consistent ly, and calculate various physical quantities. In 
Fig. [TJ we plot Tbkt with respect to the SO coupling 
strength aKp^ the Zeeman field ft, and the binding en- 
ergy E\). We find that, quite surprisingly, Tbkt decreases 
with increasing a, although the superfiuid order parame- 
ter A is enhanced [28j . The unexpected decrease of Tbkt 
has not been found in previous literature [41[ and can be 
understood as follows. For ft = and aKp <C 1, the 
superfiuid density 
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h c = y//i 2 + A 2 [1, Around ft c , the minimum quasi- 
particle energy gap occurs at k = (i.e., E g = E-^=o), 
which first closes and then reopens across ft c , allowing 
the system to change its topological order from a reg- 
ular 5-wave superfiuid to a topological superfiuid where 
MFs exist in vortex cores [8]. At finite temperature, the 
zero temperature topological phase transition becomes a 
phase crossover. In Fig. [2^, we plot the line E g = T for 
a finite ft, which is obtained by solving Eqs. (|2j), (J6j) and 
Eg = T self-consistently. The gap closes at a critical 
where ft c = vV + A 2 = ft. When E h < (>) E c h , ft > (<) 
ft c , and the region below the line E g =T corresponds to 
the topological (non-topological) superfiuid. Above the 
line Eg = T, the temperature is larger than the quasi- 
particle energy gap and the thermal excitations destroy 
the topological superfiuid. We emphasize that at finite 
temperature there is no sharp phase transition, but only 
phase crossover between topological and non-topological 
superfluids. The line E g = T is plotted only for guiding 
purpose and there is no sharp boundary between different 
phases. 

The solid and dash dotted lines in Fig. [2^ correspond 
to Tbkt and T Vor tex, respectively. We see T Vor tex quickly 
approaches a constant Tyortex = 3Ef/40tt with increas- 
ing binding energy. In the pseudogap region (T > Tbkt), 
free vortices may exist, but they are not suitable for the 
observation of MFs due to the lack of phase coherence. 
While in the vortex lattice region (T < Tyortex), the zero 
energy modes may break into two normal states with en- 
ergy splitting ex Ae~ R ^ 39], UK because of the large tun- 
neling between neighboring vortices in the lattice, where 
£ is the coherence length of the superfiuid, and R is the 
intervortex distance. In this region, strong disorder in the 
tunneling may also lead to Majorana metals [49]. Clearly, 
the required temperature for observing MFs should be 
Tyortex < T < Tbkt, instead of the intuitive zero tem- 
perature 31, 32[ . In Fig. [2^, the shadow regime between 



, (7) 



Tbkt and Tyortex gives the possible parameter range for 



decreases with increasing a. Here Eq = y^k + A 2 . 
Physically, the increased density of states near the Fermi 
surface dominates at small a, hence enhances the phase 
fluctuation and reduces the superfiuid density. Note that 
the next leading term in Q is ~ a 6 for T = and ~ a A 
for T 7^ 0, therefore the above analytical result is still 
valid even for aKp ~ Ep. The Zeeman field is detri- 
mental to the superfiuid order parameters, thus further 
reduces the superfiuid density and the BKT temperature, 
as shown in Fig. [TJd and Fig. [T]i. Perturbation theory 
near ft <C 1 shows that the change of Tbkt is ~ ft 2 , with 
the coefficient depending strongly on a. Similar features 
are also found for Tyortex- 

In the presence of both SO coupling and Zeeman 
field, a topological superfiuid can emerge from regular 
s- wave interaction when ft is larger than a critical value 
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FIG. 2: (Color online). Parameter region for topological su- 
perfluids. (a) QC: quantum critical region; TS: topological 
superfluids; NTS: non-topological superfluids. The shadow 
region is the possible parameter region for observing MFs. 
ft = 0.6Ef, aKF = 1.5Ef- (b) Phase boundary between TS 
and NTS at T = 0. The arrows mark the required minimum 
Zeeman fields. 
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51] 



the experimental observation of MFs which exist in the 
cores of the naturally present V-AV pairs in this region. 
Note that the intervortex distance R is still essential for 
the observation of MFs in this region. Such distance 
may be estimated using an analogy between V-AV pairs 
and the 2D Coulomb gases. It has been shown [54] that 
the mean-square radius (R 2 ) ~ £ 2 Zpj-l > wnicn diverges 
at T = T B kt (see Eq. (@J)) and the V-AV pair breaks 
into free vortices. Therefore there should exist a finite 
temperature regime below Tbkt where R ^> £ and the 
splitting of the zero energy Major ana states is vanish- 
ingly small. Note here that our theory can only capture 
the average behavior of V-AV pairs and a more delicate 
theory is still needed to further understand the detailed 
structure of V-AV pairs. 

In Fig. EJd, we plot the parameter region for topologi- 
cal superfluids with respect to the Zeeman field and the 
binding energy at the zero temperature. Generally, the 
required critical for topological superfluids increases 
when h increases. Note that here the phase boundary is 
determined by h c = \J \i 2 + A 2 at T = 0, but does not 
shift much even at finite temperature. 

Because MFs can only be observed at finite temper- 
ature, there exists a nonzero probability ~ exp (—rj) 
for thermal excitations to non-topological excited states 
in the vortex core, where the ratio 77 = e m /T, e 
A 2 /(2E F ) is the minimum energy gap (minigap) [5 
in the vortex core that protects zero energy MFs. For 
the observation of MFs, it is crucially important to have 
77 > 1, in addition to the requirement of the temperature 
^Vortex < T < Tbkt (i.e., 77 > t/bkt = e m /T B KT and 
V < Wortex = e m /T Vo rtex). Wnen E b > E F , we have 
?7bkt = SEb/Ep and r/vortex = 407ri^/ 'E F , which means 
77 ^> 1 for a large E b . Remarkably, 77 is also dramati- 
cally enhanced by the SO coupling in the BCS side. In 
Fig. [3j we plot 77BKT and //Vortex with respect to E b for 
parameters h = 0.8E F and aKp = 1.6E F . The vertical 
line at E b c± 0.33Ep is the boundary between topological 
and non-topological superfluids. There is a broad region 
in the BCS side with 7?bkt > 1 even at the highest tem- 
perature Tbkt • The region can be much larger when the 
temperature is further lowered to T = Tv or tex- The solid 
circle represents the possible temperature T = 0.05E F 
that ma y b e accessible in experiments in the near future 
[l2l [H, HH, yielding 77 ~ 5.6 and the thermal excita- 
tion probability less than 0.4%. Such a small probabil- 
ity clearly demonstrates that MFs can actually be ob- 
served with realistic temperature in experiments. Note 
that in the presence of SO coupling e m may be larger 
than A 2 /(2E F ) used for our estimate [51[, therefore 77 
could be even larger, which further suppresses thermal 
excitations. 

We illustrate how to find suitable parameters, in par- 
ticular Eb, for the observation of MFs. Because topolog- 
ical superfluids only exist in the region h > ^/ fi 2 + A 2 , 
we need |A| 2 ^> |/i| 2 to obtain a small h and a large 
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FIG. 3: (Color online) Plot of 77 = s m /T as a function of Eb. 
h = 0.8E F and aK F = 1.6E F . Dashed-Dotted line: El = 
0.33E F as the corresponding boundary between topological 
superfluids (TS) and non-topological superfluids (NTS) (see 
the Eg —T line in Fig. [2^, for the determination of the phase 
boundary). Dotted line: 77 = 1. Filled circle corresponds to 77 
at T = 0.05Ef. The shadow region is the possible parameter 
region for observing MFs. 



A. Clearly the BEC limit with large E^ does not work 
because \fi\ - \E F - E b /2\ > A - ^/2E b E F . While in 
the BCS side, |A| 2 ^> \ji\ 2 can be achieved by choosing 
suitable SO coupling, Zeeman field and binding energy 
[28|. To obtain a large quasiparticle minigap in the vor- 
tex core (thus a large 77) for MFs, the binding energy 
should be set to be close to the critical E% (see Fig. [3]). 
At the same time, the bulk quasiparticle gap should also 
be chosen to be much larger than the temperature. Note 
that in experiments the bulk quasiparticle gap for the 
topological superfluid can be detected using the recently 
experimentally demonstrated momentum-resolved pho- 
toemission spectroscopy [55|. 

Finally we briefly compare the cold atomic gases [9] 
with the semiconductor-superconductor nanostructure 
where certain signature of MFs has been observed in ex- 
periments [56-59]. Both systems share the same ingre- 
dients for MFs: SO coupling, Zeeman fields and s-wave 
pairing, besides the s-wave pairing is through intrinsic s- 
wave interaction in cold atoms, while externally induced 
in the nanostructure. Because of high controllability and 
free of disorder (lacked in the corresponding solid state 
systems), the topological cold atomic superfluids provide 
an ideal and promising platform for observing MFs and 
the associated non-Abelian statistics, which are both fun- 
damentally and technologically important. 

In summary, we study the BKT transition in 2D SO 
coupled Fermi superfluids and find the unexpected de- 
crease of the BKT transition and vortex lattice melting 
temperatures with increasing SO coupling. We charac- 
terize the finite temperature phase diagram for the ex- 
perimental observation of MFs in this system. Our work 
not only provides the basis for future study of rich and 
exotic 2D SO coupled Fermi superfluid physics, but also 
yields realistic parameter regions for experimentally real- 
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izing nontrivial topological superfluid states from stable 
cold atom s- wave superfluids. 
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